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Abstract 

In this note, we discuss the quantum version of the melting crystal corner in one, two, 
and three dimensions, generalizing the treatment for the quantum dimer model. Using a 
mapping to spin chains we find that the two-dimensional case (growth of random partitions) 
is integrable and leads directly to the Hamiltonian of the Heisenberg xxz ferromagnet. The 
three-dimensional case of the melting crystal corner is described in terms of a system of 
coupled xxz spin chains. We give a conjecture for its mass gap and analyze the system 
numerically. 
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1 Introduction 



In this note, we study the quantization of the statistical system of a melting crystal in one, 
two, and three dimensions. In two and three dimensions, the melting crystal configurations 
correspond to the diagrams of random partitions, resp. plane partitions. Random partitions 
appear in many branches of mathematics and physics, such as Gromov-Witten theory, ran- 
dom matrix theory, Seiberg-Witten theory and others, see e.g. [1] for an overview. Also in 
string theory, random partitions are of importance. In (21 13|, the partition function of the 
melting crystal corner given by the MacMahon function was shown to equal the partition 
function of the topological string A-model with target space X = C 3 . 

In the spirit of stochastic quantization, we construct a quantum Hamiltonian whose 
ground state corresponds to the classical steady state distribution. The treatment is reminis- 
cent of the work of Rokhsar and Kivelson on quantum dimer models |4] (the perfect match- 
ings of the dimer model on the hexagonal lattice being in one-to-one correspondence with 
the configurations of the 3d crystal) and Henley's study of the classical and quantum dynam- 
ics for systems satisfying the detailed balance condition [5[. While the quantum dimer model 
only keeps track of the total number of configurations, the number of boxes of each crystal 
melting configuration is taken into account in our system. In the partition function, each 
configuration is weighted by ^ #boxes . The quantum dimer model at the Rokhsar-Kivelson 
point corresponds therefore to the limit q — > 1 of the melting quantum crystal. 

The problem of the two-dimensional quantum crystal can be reformulated in terms of 
a one dimensional fermionic system. The system translates to an asymmetric, hard-core 
diffusion problem on a half-full infinite chain. In fact, the quantum Hamiltonian we derive 
is easily transformed into the Hamiltonian of the XXZ ferromagnetic Heisenberg spin chain. 
The ground state of the partition problem corresponds to a particular sector (the half-filled) 
of the kink ground state of the infinite chain. Using the technology of basic hypergeometric 
series, we derive the limit shape 0, which corresponds to the integral of the magnetization 
in the limit q — > 1. Also the general expression for n-point correlation functions is derived. 

Based on our findings in the two-dimensional case, we can express the 3d system as a 
system of interlaced spin chains. The quantum crystal corner is a rare example of a three- 
dimensional system which is nonetheless fairly tractable. This is of course due to its very 
constrained nature. Expressed in terms of spin chains, it becomes manifest that its dynamics 
is not even fully two-dimensional. 

Based on numerical evidence, we conjecture the mass gap of the 3d quantum crystal to 
be the same as in the Id and 2d systems. More generally, all these systems can be interpreted 
as directed random walks on the graph of configurations, where going forward and going 
back has different weights. Depending on the dimensionality of the problem, the graph of 
configurations changes, but it always remains a directed acyclic graph. This more general 
framework unifies the three problems considered here and might also account for common 
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properties, such as the value of the mass gap. 

The plan of this note is as follows. In Section |2j we explain the general procedure of 
deriving the quantum Hamiltonian for a melting crystal. As a warm-up we specialize in 
Sec. [3] to the simplest possible problem, namely Id crystal melting. The two-dimensional 
system (i.e. random partitions) is discussed in Section]!} where the Hamiltonian is expressed 



in terms of fermionic operators (Sec. |4.1|| and mapped to the XXZ spin chain (Sec. 4.2 1. The 



1-point, 2-point, and n-point correlation functions are derived in Sections |4.3.2| to 4.3.4 



while Sec. 4.4 gives some numeric results for higher eigenstates. In Section |5j the three- 



dimensional system is discussed and numeric results for the higher eigenstates of the 3d 



case are presented. The conjecture for the mass gap is stated in Sec. 5.3 



In Section |6j we end with concluding remarks and an outlook on the open problems. 

2 The quantum Hamiltonian 

In this section, we will explain the general framework of quantization we will be using. The 
quantum dimer model of Rokhsar and Kivelson [4] is an example thereof. In fact, it is a limit 
of the more general problem of quantum crystal melting. 

2.1 The general framework 

Consider a classical system with a discrete set of configurations { C, }, such as a lattice model. 
Let H be the energy functional that associates to each configuration its energy 7?(Q). If we 
consider the canonical ensemble, the probability of finding the configuration Q is given by 
its Boltzmann weight 

KQ) = — z~ / (2-D 

where /3 is the inverse temperature and Z is the partition function 

Z = £V/W) . (2.2) 

i 

We want to quantize this system. To do so we introduce the separable Hilbert space 
spanned by the orthonormal vectors |Q) which are in one-to-one correspondence with the 
classical configurations. Each vector in \ip) S J4? defines a probability measure on the Q as 
follows: 

mfiQ)= \imt, (2.3) 

where the scalar product is understood in I 2 . 
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Consider the state |Yo) defined by 

|Y ) = X>-^< )/2 |Q> • (2.4) 
i 

It is immediate to verify that the measure associated to |Yo) is precisely the Boltzmann dis- 
tribution for the canonical ensemble, up to the normalization fixed by 

<Y |Y )=Z. (2.5) 

To be more explicit, consider a physical quantity Q that can be measured on the classical 
configurations { Q }. On the canonical ensemble, the statistical average of Q is given by 

(Q)ci = ^EQ( C ')^ (Q) - (2-6) 

i 

To such a quantity we can associate an operator Q in the Hilbert space J4? which is diagonal 
in the basis given by the |Q), 

0\d) = Q(Q)|Q) . (2.7) 

For any such Q, the statistical average thus coincides with the expectation value of Q on the 
state |Yo) (which is equivalent to the above statement of the measure corresponding to the 
Boltzmann distribution): 

<X^> 1 £ [e-WWm'iQM (C ; .|Q)1 = (Q) d . (2.8) 
(t |t ; z ty 



Thus the choice of normalization in Eq. (2.5 1 is meaningful. 



Let now H be a Hamiltonian operator in the Hilbert space that admits |Yo) as its 
unique zero energy ground state, 

W|Y )=0. (2.9) 

We call Ti. the quantum Hamiltonian corresponding to 7i. Note in particular that f> is now 
just a parameter in Ji. The quantum partition function which describes the quantum behaviour 
of the system contains on the other hand a new coupling constant h and is defined as 

2£ = Tr,^(e- n/h ) . (2.10) 

For a given operator Q the expectation value is defined as the following path integral: 

(Q) = ^^MQe- n/h ). (2.11) 
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The classical limit is obtained for fi — > and using a saddle point approximation one veri- 
fies that the main contribution to the expectation value comes from the ground state, hence 
reproducing the classical result: 



(Q) = TrMQe- H/h ) — > = (Q)ci • (2-12) 



A geometric interpretation for the operator H can be found in terms of the classical con- 
figuration space. Consider the graph obtained by associating each of the C, to a node and 
linking two nodes with an edge weighted by Z;y = exp [—/J (H(Q) — 7?(Cy))] . Define the 
operator A ? on this graph as 



A, 



E Wi q (n(Ci)-H(c k ))/2 iti=j 

— 1 if i and j are connected (2.13) 

otherwise. 



where q = e $ . The notation (ft:); denotes the sum running over all configurations Q adja- 
cent to Q. One can easily verify that 

A,|T >=0. (2.14) 

Since this is the property we were looking for in the quantum Hamiltonian, we take A q = H. 
There is some arbitrariness in the choice of A q ; in fact one could have added longer range 
interactions and still preserved the annihilation of |Yo). 

The A q operator can be interpreted either as a ^-deformed Laplacian on the state graph 
or as the standard Laplacian plus a potential term: 

n/(u) = -j e \/(v) - /(«)] + / E [ 9 w«)-«w)/ 2 - il /(«) = -/ (a + v„(u)) m, 

{V)u (V) u 

(2.15) 

where / is the quantum energy scale (not to be confused with the classical coupling /3 of the 
classical Hamiltonian H), and f(u) is a function defined on the graph. The choice of the 
minus sign in front of the Hamiltonian is needed to ensure that |Yo) is the state of minimal 
energy. Note that in the jS — > limit the potential term V q is proportional to the second 
variation of the classical Hamiltonian. 

By following the above quantization description, the dynamics and a time direction 
are added to the initial classical problem, since the states in the Hilbert space satisfy the 
Schrodinger equation 

ift^|T(0>=H|T(0>. (2.16) 
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This dynamics corresponds to a random walk in the space of configurations where we intro- 
duce a (discrete) metric depending on the classical energy. 

When this article was close to completion, we became aware of flZl which has some over- 
lap with the treatment given in this section. 

2.2 The quantum dimer model and quantum crystal melting 

The quantum version of the dimer model was introduced in 1988 [4 J . Given a dimer model 
on a bipartite graph, the perfect matchings form a basis for the Hilbert space of the quantum 
dimer model (QDM). Its (quantum) Hamiltonian takes the form 

H = -JE(PX"I + !■)(□!) + vE(PXq + |")<"|) , (2.17) 

where |D) and ■) represent the two possible dimer configurations in a plaquette. 

The first term is kinetic and runs over all pairs of perfect matchings |a)(jS| which only 
differ by one elementary dimer flip move (plaquette move), in which a matching around a 
given plaquette is turned by one position, thus exchanging edges covered by a dimer with 
uncovered edges. The second, potential, term is diagonal in the Hilbert space and counts 
the number of flippable plaquettes in configuration a. The coupling constants / and V are 
phenomenological parameters. At the Rokhsar-Kivelson (rk) point V = }, the exact ground 
state wave function is an equal-weighted superposition of all perfect matchings. 

This construction can be generalized, as shown in 0. The dynamics of any discrete 
stochastic classical model is described by the master equation 

4- M?) = E i w «?W - WpM*)) > ( 2 - 18 ) 

where P«(t) is the probability to be in configuration a at time t, and W a p is the transition 
rate to state a if the system is in state f>. It was shown in fSj (see also [7J) that every discrete 
system fulfilling the detailed balance condition 

= W^pf (2.19) 

has a RK point with the same mapping of the eigenfunctions to the classical dynamics. 

The system of a melting crystal corner can be mapped to the problem of stacking cubes 
in an empty corner of 3D space. The growth rules are the following. A cube can be added 
if three of its sides will touch either the wall or the faces of other cubes. This leads to a 
minimum energy configuration without any free-standing cubes. The partition function of 
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(a) The empty room per- 
fect matching 



(b) After one plaquette 
move 



(c) After five plaquette 
moves 



Figure 1: Empty room perfect matching and matchings after flip moves 



the melting crystal corner takes the following form: 

00 

Z= £ q N & = Yl(l-q n )-\ (2.20) 

3d partitions a n—1 

where N(a) is the number of boxes in the configuration a. and the rightmost expression 
is the so-called MacMahon function (8). The convergence of the above expression is only 
guaranteed for < q < 1, which we will assume henceforth. 

The configurations of the melting crystal corner are in one-to-one correspondence with 
the perfect matchings on an infinite hexagonal lattice. There is a unique (up to translation) 
empty room perfect matching which serves as the starting point. A plaquette move corre- 
sponds to adding or removing a box (21, see FigureJTJ Nonetheless, the QDM explained above 
does not correspond directly to the quantum system for the melting crystal corner, since the 
number of boxes does not enter anywhere. In the following, we derive the quantum Hamil- 
tonian for the crystal melting. Even though the following paragraph phrases everything in 
the language of the three-dimensional crystal melting problem, the procedure is completely 
general and applicable to any dimension. 

To describe the (classical) dynamics of the crystal melting process, we express it as a 
random walk on the state graph (where each node corresponds to a configuration and two 
nodes are joined by a directed line if the endpoint can be reached from the initial one by 
adding a box - see Fig.|2|. The dynamics is encoded in a master equation of the type ( 2.18| . 

We choose the following transition rates: 



W ajS 



if N(j8) =N(a) + l, 
y/q ifN(0) = N(a)-l, (2-21) 
otherwise. 
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Figure 2: State graph of the melting crystal corner for the first three levels 



In other words, 

W^C^M-MW)/*, (2.22) 

where is the adjacency matrix for the state graph. One can verify that the (unique) 
stationary distribution for P a is 

Pf = \ q NM , (2.23) 

where, as in Eq. ( 2.20| l, 

Z = Y j q N{ * ) - (2-24) 

a 

This is equivalent to saying that the system is in a Boltzmann distribution with classical 
energy 

H{ol) =N(a). (2.25) 

In fact, one can easily verify not only that P« is an equilibrium state, but that it also satisfies 
the detailed balance condition ( 2.19) . It is customary to define the exit rate from state a as 

W aa = - 7a = - £ = - £ C a ^- N ^ /2 = deg» - ^ deg-(«) , 

(2.26) 

where deg^(u) is the number of configurations connected to a containing one block more 
(less) than a. Note that in this way, the evolution can be described by the vector equation 

^P(0=WP(t), (2.27) 
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and the stationarity condition becomes 

WP (0) = _ 



(2.28) 



Instead of using the matrix W, one can define the symmetrized version 

W a p = (Pi 0) ) ~ m W a/5 (pPy U1 (2.29) 
(no summation implied). Explicitly one finds 

W a « = G»« if ot ^ B 

y P P (2.30) 

W M = W aa . 

Given the characteristic equation W u p<f>^ = one can verify that the ground state is 

now represented by the vector ffi = \J . 

Let us now switch to the | Denotation used in the QDM. The ket ■) here represents a 
cube in the configuration that can be removed, while |D) represents a place where a cube 
can be added. Now, W aa can be written as 



W aa = (a | 



£^|n><q + -i=|B)<i 



\cc), (2.31) 



where the sum runs over all places where a cube can be added or removed in the configura- 
tion. Similarly, the non-diagonal terms can be written as 

w^ = <«iE !■)(□! + !□)(■!] 1/3). (2.32) 

It is therefore natural to define the quantum Hamiltonian by 

n = -/El")<q + !□)(■! + !□>(□! + -J= ■ ■ , (2.33) 

acting on the Hilbert space generated by the orthonormal basis of the configurations \oc). 
Notice that at the RK point V/J = 1, this coincides with W: 

H = -/W. (2.34) 

Being proportional, 7i and W share the same eigenvectors with proportional eigenvalues: 

n^4>^ = -J\4>i X) . (2.35) 
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In particular, the ground state Eq = corresponds to 



tfW^^^ 72 ' (2-36) 
which has to be interpreted now in the sense of a quantum superposition: 

(ground) =Y j q N{a),1 W) , (2.37) 

where |#) is the state corresponding to the perfect matching a. Note that we changed the 
normalization for the wave function as in Eq. \2.5) , so that 

a 



ground|ground) = J^^^ = z ■ ( 2 -38) 



As as consequence of the detailed balance condition in Eq. ( 2.191 , the ground state is 
frustration free, i.e. each local interaction is minimized separately. 



3 The one dimensional problem 

As a warm-up we briefly consider the one-dimensional analog to the crystal melting prob- 
lem, i.e. a system where a set of blocks is added or removed along a line with integer lattice 
points starting from n = 0. The states are labelled by an integer number \n) and the quantum 
Hamiltonian acts on them as follows (note that we placed ourselves at the V = } point): 



« |0> = -/ [|1> - f 1 ' 2 |0>] , 

U \n) = -} [\n + 1) + \n - 1) - (q 1/2 + q' 1 ' 2 ) \n)] , V n > 



2 

m 



(3.1) 



As shown in the previous section, the ground state is 

CO CO -1 

|Yo> = E V" /2 \ n ) > Z = Wo) = £ q» = :— • (3.2) 

(3=0 n = 1 " 

For n > 0, the Hamiltonian can be written in terms of the Laplacian on the line as 

H = -jA + m 2 , (3.3) 

where 

A|n) = \n + l) + \n-l)-2\n) , (3.4) 



j(q 1/2 + q- 1/2 -2) . (3.5) 
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Hence we can immediately read off the value for the mass gap. In fact the system is easily 
solvable and one finds the eigenvectors 

|Y(*)) = E (sin(kn) + q 1/2 sm(k(n + l))) \n) , fc^O, (3.6) 

J! = 

with energies 

H \Y(k)) = } [q in + q- in - 2cosfc) |T(Jfc)> . (3.7) 

This set of eigenstates does not contain the first excited state |Yi) which would correspond 
to k = 0. It satisfies 

H |Yi) = m 1 |Yi) . (3.8) 

This means that |Yi) is harmonic: 

A|Yi)=0. (3.9) 

In one dimension, this is easily solved and using the equation for |0) as initial condition one 
finds 

l Y i> = E (-n + q 1/2 (n + l)) \n) . (3.10) 

n=0 

4 The two-dimensional problem: quantum partitions 

After having discussed the simplest version of the problem in the last section, we now turn 
to the study of the growth of a two-dimensional crystal. When using the picture of squares 
being added to an empty corner of the plane, the growth rules are the following: a square 
can only be added when two of its sides will touch either the wall or the side of another 
square. It is easy to see that in this case the allowed crystal configurations are in one-to-one 
correspondence with random partitions. A partition of n is a non-increasing finite sequence 
of positive integers whose sum is n. We can visually represent a partition by its (Ferrers) 
diagram, for which we adhere to the Russian tradition, see Figure |3J where the corner of the 
plane is rotated by 135 degrees. The theory of partitions is at the same time an old subject 
of mathematics with many important results dating back to Euler, and one which has seen 
major progress in the last fifty years (see e.g. [9| for an overview). 

Let { C, } be the set of all integer partitions (or equivalently Young diagrams). Since we 
are interested in the size of a partition, it is natural to define the "classical" Hamiltonian 

7?[Q] = # of boxes of the partition Q. (4.1) 
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The classical partition function for 2d partitions is given by 



Z 2d = £ q H[Q] = Yl(l- q")-\ (4.2) 
Q n=l 

Using the procedure detailed in Sec.|2]we define the graph of configurations whose nodes 
are in one-to-one correspondence with the Q and where a line is drawn between two parti- 
tions if they differ by only one square. The Hamiltonian retains the form given in Eq. ( 2.33} , 



n = -/£ !■)(□! + !□)(■! + v£ vWX n l + -L ■ ■ , (4.3) 

acting on the Hilbert space generated by the orthonormal basis of the configurations |Q). 
The ket ■} now represents a square that can be removed from the partition, while (□) 
represents a place where a square can be added. The sum runs over all places where a 
square can be added or removed in the configuration. 

The 2d case of random partitions is especially interesting for us because of the existence 
of a map between 2d partitions and the Neveu-Schwarz sector of the complex fermionic 
oscillator. It enables us to express the quantum Hamiltonian ( |4.3| in terms of fermionic 
operators. 

4.1 Fermion operator formalism for 2d partitions 

In this section, we will explain the map from partitions to fermions mentioned above. The 
empty partition is depicted in Figure|3j(a). The horizontal line shows the projection to the Id 
fermion problem, where the fermions sit at half-integer positions. All the negative positions 
are filled, therefore the empty partition corresponds to the filled Fermi sea. The fermion 
operators obey the usual anti-commutation relations 

{</>«, CI = Km, n,m£Z + 1/2, (4.4) 

all other anti-commutators zero, and the annihilation operators annihilate the vacuum, 

xp n \Q)=0 VneZ + 1/2. (4.5) 

Note that the vacuum |0) corresponds to the bottom of the Fermi sea. In terms of fermions, 
the empty configuration (no squares) is represented by 

^ oo ...^*i|0) = |haIf) / (4.6) 
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(a) Empty partition 



(b) Partition with one box 



Figure 3: Empty partition and partition with one square with projection to Id fermions 

where each NW-SE diagonal line corresponds to a creator tp*. The state |half), corresponding 
to the filled Fermi sea (a picture also familiar from matrix models), obeys 

xp*_ n |half) = 0, \p n |half) =0, n > 0. (4.7) 

The partition with one square, see Figure [3] (b), corresponds to 

tp_ 3 ip\ |half) = tp*_ co ...xp*_ 3 xp* 1 \0), (4.8) 

2 2 2 2 

which corresponds to an annihilator for the SW-NE line at — \ and a creator for the NW-SE 
line at \. On the projection to the horizontal line, the fermion at position — | was annihilated 
and a fermion was created at j, i.e. the fermion at — ^ hopped over to \. Adding squares 
to the partition corresponds to fermions hopping to the right, in this sense this can be seen 
as a diffusion process in which every position can be occupied only by one particle, {i.e. it 
is exclusive). Figure [4] shows a partition with the corresponding occupied and free positions 
shown in the projection to the horizontal line. The Hamiltonian (|4.3|l becomes now 



/ V VI \ 

n 2d = -J\ 2j ^m+i^Pm + 1/41/W1 - T ^qn m (1 -n m+1 ) - -— n m+l (1 -n m ) ), 

\meZ+l/2 ' > V" J 

(4.9) 

where n m = ipl n ip m is the fermion number operator. The first two terms are the kinetic hop- 
ping terms, while the last two terms are the potential terms, counting the number of possi- 
bilities to hop right and left per configuration. The different weights for hopping right and 
left ( ^Jq and 1 / ^/q) make the system into an asymmetric diffusion process. The fermion num- 
ber rij is conserved in this system, therefore states with differing fermion numbers belong to 
different superselection sectors S„ f of the Hamiltonian. 

Let us now go to the RK point V = }. The unique ground state, satisfying H | ground) = 0, 
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lt>_ 15 lb_9 llj_5 tp*l 1p% 1p*19 



Figure 4: Two-dimensional partition corresponding to the state tpl? ip_ 15 9 tpl xp_ s |half) . 

2 2 2 2 2 2 



is written as in Eq. ( |2.37| : 



| ground) = lim 



E 



ai<U2<---<a l i 
bi<b 2 <--<b d 



z=l 



|half) , 



(4.10) 



where we notice that the number of squares in a given partition }i is given by 



d 

I 

!=1 



(4.11) 



4.2 The Heisenberg xxz spin chain 



The Hamiltonian in Eq. (4.9 1 can be recast into a different form, which is more familiar in 



statistical physics. If we express the creation and annihilation operators in terms of the usual 
Pauli matrices a^, k = 1,2,3 at position m, Eq. (4.9 1 becomes^ 



H 



2d 



/ 



E 



meZ+l/2 



1 



(4.12) 



Going to the RK point V/J = 1, this is precisely the Hamiltonian for the Heisenberg xxz 
spin chain |10| with anisotropy parameter 



A 



(4.13) 



1 Note that our notation differs form the one generally used in condensed matter literature, where not Jq, but 
q' is used which is related to q by q' = Jq. 
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Note that the term proportional to (cr^ — c^+i) * s * ne on ly term in Eq. (4.12 not invariant 
under the exchange q — > 1/q. It is a boundary contribution vanishing for the case of the 
infinite chain. In fact one could have noticed this symmetry before since for any given 2d 
partition one square more can be added than removed. It follows that the asymmetry be- 
tween q and 1 / q is only proportional to a trivial constant term. This implies also that we can 
restrict our attention in this case to V = J without loss of generality. 

Since for the partitions, we are interested in the case < q < 1, we have 1 < A < oo, 
which means that we are in the ferromagnetic regime. The case A = 1 corresponds to the 
isotropic Heisenberg or XXX model. For A = oo, this is the one-dimensional Ising model. 
The xxz model is one of the best studied systems in statistical mechanics and can be solved 
using the Bethe ansatz [11 j. It was shown to correspond to the six-vertex or ice-type model. 
That the finite XXZ chain admits the quantum group U ^[SU(2)} as a symmetry was first 
pointed out by Pasquier and Saleur |T2). A large part of the literature is devoted to the 
anti-ferromagnetic case or the case — 1 < A < 1, but a fair amount of results exists on the 
ferromagnet, which is collected in Ifl3l . 

We are specifically interested in the infinite volume chain. It admits four types of zero 
energy ground states. Two are translation invariant and correspond to all spins up and all 
spins down. The other two are the kink and the anti-kink fl4lll5] . This list of ground states 
was shown to be complete. These ground states are furthermore frustration free, i.e. they 
minimize each of the next-neighbour interactions separately. Over each of these ground 
states, a spectral gap exists with value IH6H 

l=-l{q 1/2 + q- in -2) ■ (4.14) 

It was furthermore shown in fFT} that droplet excitations exist over all four ground states. In 
our normalization, their energies are 

£(") = -'(^-1 ,/2 )^- C«5> 

where n is the number of spins that deviate from the ground state and g = — log q. For 



n = 0, £(0) = 0, while for n = 1, the mass gap Eq. (4.14 is reproduced. These energy levels 



are symmetric under the exchange q — > 1 / q as shown in Fig. |5j 



4.3 Correlation functions for the kink states 

Since the partitions correspond to the case of the kink centered in 0, we concentrate on the 
kink states. The kinks form an infinite family of ground states interpolating between spin 
up at —oo and spin down at +oo. These kinks are the lowest energy eigenstates for the 
superselection sectors with fixed particle number into which the Hilbert space separates 
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Figure 5: Energy levels for the droplet states as a function of f> for different values of the spin 
deviation as in Eq. ( 4.15| |. 



(note that [H, £ m cr^ ] = 0, which means that the action of the Hamiltonian does not change 
the particle number). 

4.3.1 Grand canonical ground state 

A useful representation for the kink state can be obtained in terms of the grand canonical 
ground state 

OO 00 y \ 

Y(z) = Ft f 1 + z ^ n+1/2 r n ) |0) = ( J.+1/2 . (4-16) 

which represents a kink centered in xq = — log |z| / log q. Note that we shifted the spin chain 
such that the leftmost spin is at n = 1/2. The kink is exponentially localized |18| . In fact if 
one considers the magnetization profile w%(z) {x) = (Y(z) |af |Y(z)), one finds that 

2 - m Y(z) 0) < ^ M iix<x Q 
\ + m Y(2) (x) < ^!' T - x of if x > x . 

For q — > 0, the system becomes an Ising model, where the kink is localized in a single 
point, i.e. the magnetization profile becomes a step function centered in x = Xq. For q = 1, 
where the symmetry is the usual SU(2), the system is the XXX spin chain which only admits 
translation-invariant ground states (one could say that the kink is completely delocalized). 



The grand-canonical kink (4.16 1 can be seen as a generating function for the kinks in the 
N-particle sectors: 

|Y(z)) = lim £ |Y N )z N . (4.18) 

In particular, our ground state corresponds to the half-filled infinite spin chain, i.e. the kink 
centered in the middle of the chain and is therefore precisely \*¥l/2)- 
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Our strategy to compute correlation functions is the following: given an operator 6 that 
does not depend on z and leaves the fermion number invariant, one can expand the correla- 
tor on the grand-canonical kink as 

oo 

<Y(z)|€?|Y(z)) = £ z 2N <Y N |^|Y N ) , (4.19) 

N=0 

where we used the fact that the |Y N ) are orthogonal since they describe states with different 
particle numbers. It follows that the correlation function on the N-particle kink is given by 

(Y M |^Y N ) ^ 1 f dz (Y(z)HY(z)) 
{ )n ~ <Y n |Y n ) - <Y N |Y N > f 2m z^+i " ' ^2u) 

This is particularly useful when does not mix contributions from different points in which 
case the correlator on the grand-canonical kink can be found easily. 
As a first example consider & = 1, i.e. the norm 

<Y(z)|Y(z)) = n (l +^V +1/2 ) = (~wq 1/2 ;qU, (4.21) 

where we used the notation of the ^-shifted factorials and w = z 2 . Using the ^-binomial 
theorem, 

oo n N 2 /2 

(-wq l/1 ;q) = £ -L_ H ; N / (4.22) 
N=0 (<MJn 



we immediately find that the norm of an N-particle kink is 
where the factor q N2 ^ 2 takes the shift of the spin chain into account. 



4.3.2 One-point function 

In general, other quantities of interest such as n-point correlation functions can be evaluated 
for Y(z) and then expanded in a series of w to get the contribution of the N-particle kink. 
As an example, let ]i(x) denote the profile of a partition ]d, (i.e. a piecewise C°° function). 
The map to fermionic states is such that if }i'(x) = — 1, there is a spin up at position x in the 
chain, while if fi'(x) = 1, there is a spin down at x . Consider the operator of acting on the 
point x, 

^IT> = ^IT>, 4\l) = -\\D- (4-24) 
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If we define the magnetization for a state |0) by 
it is clear that 

1 d 

m^{x) = --—^{x), (4.26) 

where }iq>(x) is the shape of the partition associated to |<E>). The evaluation of the one-point 
function for the generating function of the kinks yields 

^)i^w> = 2 ?r + ^) n('+^ 1/2 ) • «*> 

By developing in a series (see App. [Ajl one can show that the contribution of the N-particle 
kink is given by 

«*(*) = (^)n = ^f^f =~U E(^ N ; *U KX+1) ■ (4-28) 



2 k= 



Shifting the center and sending N — > oo one finds the exact expression for the shape of the 
kink: 

1 00 

m4*) = ~ + £(-l)V 2 V"- (4-29) 
z k=0 

To study the q — > l - limit one can introduce a ^-difference equation satisfied by WooOO, see 
App.|A] It follows that the magnetization profile is 

w(ac) = -^tanh(|). (4.30) 

The shape of the partition is hence given by 

H^x) = -2 J dxm^x) = 21og(2cosh(f )) (4.31) 

or, in a more symmetric form (take ji,p{x) =: y), 

e (x-y)/2 + e -(x+y)/2 = 1/ ( 432 ) 

see Fig. [6] This reproduces the limit shape for 2d partitions given in J6). 
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(a) Magnetization 



(b) Limit shape 



Figure 6: Magnetization (a) and limit shape (b) for the ground state of the two-dimensional 
partition. 

4.3.3 Two-point function 

Given the factorized form of Y(z) in Eq. ( 4.16| l, one can easily evaluate two-point functions. 
In particular, 

Let q x ' = Z,{. We have the expansion 

We can now repeat the same construction as for the one-point function and find the contri- 
bution of the N-particle kink: 

KOn- (Yn | Yn) "4 + 2(l-^)^ I 9 9 J" (435) 

In the N — > oo limit, this becomes 

<<£<*>- = J + 2^S? £(-!)' (?'*■ ■ 

To take the ^ — > 1~ limit, we define /(x) as in Eq. ( |A.15| and, after the rescaling X{ — > 
-xj log fa), 

3 3\ 1 1 — x i \ ( e Xl e * 2 



W- ~ 4 " 2 COth (^— ) V 1 + ^ 1 + 



i tanh(|) tanh(|) = , (4.37) 
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which shows that the two-point correlation function factorizes. The same holds true for the 
n-point functions as we show in the following section. 



4.3.4 The n-point function 

The calculation of the n-point function can be simplified by introducing a recursive proce- 
dure to obtain the (n + 1) -point function from the n-point function. 
The n-point function for the grand-canonical ensemble is written as 

(Y( z )|4...4|Y(z)) = in^^(- w ^2.^ )oo = / («) (Ci Cn ) { - aq m ;qU 

(4.38) 

where w = z 2 and = q Xi+l ' 2 . It is convenient to develop it in a series in w. 

CO 

Q = £>*A£%i Cn). (4-39) 

Passing to the (n + 1) -point function just amounts to multiplying by an extra factor of 

1-Z"Ch+1 . 

~r(n+l), r ~ ~ s _ 1 — r(»)/-r r \ 

f/w Ifcl/- • -/fen/fen+lj — 1 , ~7 Jw {{,!,- ■ -,{,n) ■ 

l-rw C, n +l 



(4.40) 



Developing in a series, 



one finds 



where 



], ^" +1 =-l + 2E(-lf^ +1 , (4.41) 



/;1T +1) (Ci In, Cn+l) = E (Ci , (4-42) 

fc=0 



A^Cfc C« + i) = -^(Ci U + (-l) fc ^ + iE(-l) ; Aj n) (Ci WCi^ 

7=0 

(4.43) 

which gives recursively the in + 1) -point function for the canonical ensemble a[" +1 ^ in 
terms of the n-point function. The recursion relation can be solved using as initial condi- 
tion the one-point function found in App.[A| 



otherwise. 
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and yields 



1 



a! " ,(?1 M=\'\t(n^)^ « (4 ' 45) 

where the indices i and / are understood in Z n . 

Now we can again use the same procedure outlined in App.|A]and find that on the N- 
particle kink, 

K-Os-^ + ^t DC; ih r t (rj ■ («e 

and in the N — > oo limit this becomes: 

1 1 oo n ( n-l n Xi i „X; + i \ 

(4 • • • <£>. = + 2^1 E ("I) 7 9® E (EI j ■ (4-47) 

In the q — > 1 _ limit, this is 

/ 3 3 \ 1 1 V- /tT A^i-*^/^ ^ 

« ■ ■ ■ = " e c ° th ( — 2~^) J = 



We see that the system decorrelates completely. 



7^fltanh(|)=n(4) oo . (4.48) 



4.4 Height function and numerics for higher eigenstates 

We would like to form an idea about the excited eigenstates of the 2d quantum crystal. Since 
at this point, we do not have an exact expression, we resort to numerical calculations for 
finite systems containing partitions up to a certain number. 

The height function is defined on the endpoints of the squares making up the partitions. 
For a system allowing partitions fitting in a box of length No, the height function is repre- 
sented on a Id integer lattice of length 2No + 1. We can visually represent any eigenstate 
by plotting its average height. So we first have to express each individual partition in terms 
of a height function and for a given eigenstate sum all these heights with the probability 
coefficients of this eigenstate. More precisely, after the quantization of the system the height 
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function becomes an observable and as such it has an expectation value for a given state |Y), 
(Mn)) Y = (Y '^' Y) = -^E^W |<Y|*>| 2 = ^(nKW, (4.49) 

\ I / ||i||a a 

where the sum extends over the partitions \a). 

We consider the contribution of each partition to the height function minus the profile 
of the empty corner, which has the form { . . ., 2, 1, 0, 1, 2, . . . }. A partition of N is given 
by N positive integers (ao, a\, . . . , flw-i) with X^q 1 a i = Each a, gives rise to a vector S aj 
with all entries zero except for a sequence of a, twos from position i — a, + 1 to i. The sum of 
these N contributions gives the profile of the partition. The height function of a partition y. 
at position m is thus given by 



h v (m) = h (m)+2 £ X[i-a i+ i,i\(™) = h (m) + S p (m) , (4.50) 



;'=0 

where ho is the height function of the empty corner and Xk,b] (w) is the characteristic function 
of the interval [a, b]. To make an example, the partition (3,2,2) gives rise to the profile 

&(3A2) = {3,4,3,4,5,4,3} (4.51) 

consisting of the profile of the empty partition ft (o,o,o) = {3/2, 1, 0, 1, 2, 3 } and 

<*(3,2,2) = {0,2,2,4,4,2,0} 
= {0,2,2,2,0,0,0} 

(4.52) 

+ {0,0,0,2,2,0,0} 
+ {0,0,0,0,2,2,0}. 

The average heights (i.e. the integrals of the average magnetization of the underlying 
spin chains) of the ground state, first, fourth and twelfth eigenstates of a system containing 
partitions of integers up to 20 (which corresponds to 2714 partitions) for a value of q = 0.9 
are depicted in Figure [7| 

' see 



The ground state corresponds to the curve we have calculated explicitly in Sec. 4.3.2 
Figure |6jb). The profiles of the next excited states are a bit flatter. Higher excited states such 
as the twelfth shown here can show a deviation from this pattern. 
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Figure 7: Average heights for eigenstates of a system of partitions of up to 20 (corresponding 
to a total of 2714 partitions) and q = 0.9. Ground state (blue), first excited state (yellow), 
fourth excited state (red), twelfth excited state (green). 



5 The 3d Hamiltonian 

5.1 Fermionic formalism 

The general problem of the 3d melting crystal corner has already been explained in Sec- 



tion 2.2 We now want to express the 3d quantum Hamiltonian in terms of fermionic opera- 



tors, analogously to the treatment of the 2d case in Section 4.1 

Not surprisingly, the quantum Hamiltonian for the three-dimensional partitions can be 
recast in terms of a system of Heisenberg xxz chains. To do so, we first map the three- 
dimensional partitions into a system of vicious walkers. Vicious walkers were introduced 
by Fisher in [19] as a system of N particles moving in non-intersecting random walks. The 
mapping to plane partitions is obtained as follows. Consider the finite problem in a cubic box 
of length N. The empty partition can be represented as an equilateral hexagon tessellated 
with rhombi. Choose the midpoints of the sides of these rhombi on the lower left side of 
the hexagon and join them with the corresponding midpoints of the rhombi on the upper 
right edge, following the sides of the cubes (Fig. |8pb)). One can verify that each partition 
corresponds to a realization of N non-intersecting directed paths, i.e. directed vicious walks 
on the lattice obtained by taking the midpoints of the rhombi in the tessellation of the empty 
partition. In particular, the empty partition is obtained as a set of parallel walks that first 
move to the right and then up. 

It is convenient to represent the system on a square lattice by tilting the picture (Fig.^c)). 
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One can easily see that the shape of each vicious walk corresponds to a Young diagram ob- 
tained by parallel-slicing the three-dimensional partition. The set of N paths can be further 
mapped to a two-dimensional system of fermions living on the vertices of the lattice joined 
by a path (Fig.[8]jd)). Let us now consider the quantum Hamiltonian of this system. Adding 
a cube corresponds to moving a fermion from to [i — 1,7 + 1). But this can only happen 
if the points (z — and + 1) are already occupied. Similarly removing a cube corre- 
sponds to jumping in the opposite way respecting the same conditions. The Hamiltonian 
can therefore be written as 

n 3d = -/£ W^ ( *- y+1) ni-i^y + i , (5.1) 
hi 

where 

= . + y*.y kl _ q m n .. (1 _ nkj) _ q -m nKi _ n . .) (5 2) 

is the Hamiltonian of the two-dimensional problerrj^] We find that the Hamiltonian de- 
scribes an infinite system of parallel XXZ chains, each interacting with the two neighbouring 
ones. Equivalently, the system also represents the dynamics of a heap of dimers |20j |2T| . 
Interestingly enough, the same Hamiltonian is obtained when considering a diagonal slic- 
ing as shown in Fig.[9]instead of the vicious walker description which amounts to a parallel 
slicing^} 

There is a direct mapping between the rhombus tiling making up the shape of the cubes 
(see Figure[8pb)) and the dimer model on the hexagonal graph (see also FigureJTJ. To the long 
diagonal of each rhombus we associate an edge covered by a dimer in the underlying lattice. 
In this way, a perfect matching determines the orientations of the various rhombi in the tiling 
and therefore the shape of the 3d partition. From the dimer picture on the hexagonal lattice 
we can arrive directly at the fermion picture in Figure [9jjc) by drawing a black dot for each 
horizontal line covered by a dimer and drawing a white dot for each uncovered horizontal 
line on the hexagonal graph. This also leads directly to the half-integer shift between adja- 
cent spin chains. Furthermore, the dependence on the two neighboring spin chains becomes 
manifest, since information about the horizontal lines of the neighboring cells is needed to 
reconstruct the complete perfect matching from the knowledge of the dimers on the hori- 
zontal lines only. 



2 In the study of the three-dimensional system we restrict ourselves to the RK point } = V. Note that while 
it is easy to generalize the lower-dimensional models (as already commented after Eq. 4.12} , in the case at hand 
the parameter V cannot be reabsorbed by redefining A. 

3 In this diagonal slicing description, the black dots in Fig.|9](c) are allowed to move up if not only the next 
position on the same vertical line is free, but also the next positions above on both neighbouring spin chains are 
unoccupied. 
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(a) Cubes (b) Vicious walkers (c) Vicious walkers (d) 2d fermions 

on the square lattice 



Figure 8: Vicious walkers and three-dimensional partitions. Representation of a three- 
dimensional partition as cubes (a); as a system of vicious walkers (b); as vicious walkers 
in a square lattice (c) and as a system of 2d fermions (d). 




(a) Cubes (b) Diagonal slicing (c) Lattice fermions 



Figure 9: Representation of a three dimensional partition as cubes (a); slicing into 2d parti- 
tions (b); lattice fermions (c). In the empty partition case the fermions occupy all the points 
below the diagonal lines. 
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5.2 Hexagonal lattice description 

We can make use of the correspondence to the dimer model on the hexagonal lattice and 
describe the fermions in the plane more conveniently using an overcomplete basis. The pla- 
quettes of a planar hexagonal lattice are labeled by a triplet of positive integers { n\, n 2/ H3 } 
such that the center of each hexagon has coordinates n\e\ + n 2 e 2 + "3^3. The e ; are unit 
vectors that can be chosen to be 

ei = (- cos f , - sin f ) , e 2 = (cos f , - sin § ) , and e 3 = (0 / l). (5.3) 

Note that the three unit vectors are not linearly independent and satisfy ei + e 2 + e3 = 
which implies that two triplets { n\, n-i, } and { m\, mi, 1113 } represent the same point if 
Ei = i (mi — rii) e, = 0. We thus obtain a system where each fermion can jump from the point 
x to x + e3 if both the points x — ei and x — e 2 are free (see Fig.^c)). Accordingly we rewrite 
the Hamiltonian in Eq. ( |5.1| as 

^3d = -/E W 2V X+e3 (! - «x-e a ) (1 " nx-e 2 ) • (5.4) 
x 

In this notation, the empty partition corresponds to the state 

oo 

|half)= n Ke 1+ » 2 e 2 |0). (5.5) 

A generic partition is expressed by 

l^)=n^x i + 7i e 3 ^x i |half) / (5.6) 

i 

where x, = a,-ei + jS ? e2 and a,-, f>{, and 7, are all positive. To be admissible, a set { x,-, 7, } 
has to satisfy certain interlacing conditions to make sure it defines a valid plane partition. 
These conditions are particularly simple in this notation once one realizes that 7;) 
are the coordinates of one of the corners of the topmost box in each column forming the 
three-dimensional partition. This means that for each column we have to verify that the two 
columns behind are taller, namely: 

7i < 7j if Xj = Xj + ei or x,- = Xj + e 2 . (5.7) 

The number of boxes in a given partition }i is 

W=E7«- (5-8) 
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The ground state of the quantum crystal expressed in fermionic language is given by 



(ground) 



.{ *i,7i } 1 



|half) , 



(5.9) 



where the sum runs over all the sets { X\, 7/ } satisfying the condition given in Eq. (|5.7| 



This representation permits an explicit mapping between three-dimensional partitions 
and perfect matchings on the (bipartite) hexagonal lattice. Each perfect matching PM defines 
a unique height function on the plaquettes, up to a constant. This function hpM is defined as 
follows: 

/ ±1 if eePM, 



h PM (A) - h PM (B) 



T2 if e £ PM, 



(5.10) 



where A and B are two plaquettes that share the edge e, the sign depends on the orientation^] 
The empty room height function corresponds to the unique (up to translation) empty room 
perfect matching (See Figure [l](a)) and is given by 



/io(«iei +n 2 e 2 + n 3 e 3 ) = n 1 + n 2 + n 3 - 3mm{n 1 ,n 2 ,n 3 } + k, 



(5.11) 



where A: is a constant which we set here to k = 0. 



The map from the partition \i expressed as in Eq. (5.6 1 to the height function h^(y) and 
hence a perfect matching is obtained as follows: 



7.-1 

MY) =My)+3£E S(y - Xi - ne 3 ) = h (y) + ^(y) , 
i n=0 



where 



S(x) 



1 if x = 0, 
otherwise. 



Using these normalizations, 



3 y 



(5.12) 



(5.13) 



(5.14) 



Numerical results for higher eigenstates 

As already explained in Sec. |4.4| we are interested in the average height function for the 
excited states. An analytic expression is at this point beyond our reach. We have nevertheless 
run some numeric simulations for a system containing plane partitions of integers of up to 

4 Given the bipartite nature of the hexagonal lattice, we can divide vertices in black and white. The sign 
depends on whether e is crossed from A to B with a black or a white node on the left. 
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12 (resulting in a total of 3462 three-dimensional configurations). Some of the results are 



shown in Figure 10 In particular one finds that the excited levels are more spread out than 
the ground state and we have some evidence that the first excited level has Z3 symmetry 
and is non-degenerate, while the second excited level has only Z2 symmetry and is thrice 
degenerate. 

5.3 The mass gap 



We have remarked already that the mass gaps of the 1-dimensional (see formula ( 3.4 1) and 



two-dimensional quantum crystals (see Eq. ( |4.14 |) are the same. For the 3d-quantum crystal, 



we have strong numerical evidence for the mass gap taking again at the same value, see 



Figure 11 We therefore conjecture that 



The mass gap of the quantum crystal is the same in all dimensions and equals to 

7 = -/(2-v?-i/V?). 

Why this should be true can heuristically be understood as follows. Consider the graph of 
configurations (Fig. [2]), which is a directed, acyclic graph. This graph admits a height func- 
tion, the height of a configuration being its number of boxes. The quantum crystal melting 
can be considered an asymmetric random walk on this graph, where a particle hops from 
one configuration to the next. This random walk corresponds directly to the 1-particle sector 
of an XXZ model living on the state graph. In one dimension, this corresponds directly to the 
XXZ model on the line with one particle moving around. In two dimensions, it corresponds 
to the half-filled sector of the Id xxz model and in three dimensions to the quantum crystal. 
In all cases the Hamiltonian takes the form 



n = -j 



A - [q 1/2 - 1 ) deg~(a) - (q~ 1/2 - 1 ) deg + (a) , (5.15) 



where deg~(a) and deg + (a) are the outdegree and indegree of the node a. We would like 
to argue that the mass gap is the same for the 1-particle sector of the xxz model on any 
directed, acyclic graph satisfying the condition 

^CM , 1. (5.16) 

deg {(X) deg~(a)^oo 

This is the case in the one-dimensional system in Sec. |3](deg~(a) = deg + (a) = 1) and the 
two-dimensional system of Sec. [I] (deg - (a) = deg + (a) + 1). We conjecture it to be the same 
in three dimensions for the quantum crystal. 

A similar result was already proven in [22J for the case of the graph of configurations 
being a tree. 
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(a) Ground state (b) First excited (c) Second excited (d) Second excited 

(degenerate) 





Figure 11: The mass gap as a function of ^fq. The red line corresponds toy = q ' + q ' — 2. 
The blue points result from a numerical simulation in a cube of size 3x3x4 with a total of 
4116 explored three-dimensional partitions. 
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6 Conclusions 



In this note we have studied systems of crystal melting in one, two and three dimensions. 
Generalizing the work of Rokhsar and Kivelson on the quantum dimer we obtain a quantum 
Hamiltonian whose unique ground state reproduces the classical Boltzmann distribution for 
the canonical ensemble. 

The Hamiltonian becomes tractable once one realizes that it is possible to introduce a 
mapping to spin systems. More precisely we found that the two-dimensional quantum 
crystal problem (which can also be phrased in terms of random partitions) is integrable, since 
it corresponds directly to the half-full sector of the XXZ ferromagnetic spin chain. In particu- 
lar, the random height function can be interpreted as the integral of the magnetization of the 
underlying spin system. In a similar fashion, the three-dimensional quantum crystal can be 
described in terms of an infinite system of coupled XXZ spin chains. This, in turn, allowed us 
to carry out a numerical analysis, leading to the conjecture of a mass gap whose value does 
not depend on the dimensionality of the problem. 

From this point on there is a variety of interesting questions to address, mainly concern- 
ing the three dimensional system. A rigorous proof for the value of the mass gap is required. 
We would also like to obtain a better understanding of the full spectrum and explicit expres- 
sions for the higher-order eigenvectors in terms of three-dimensional partitions. Based on 
the number of conserved quantities, one is inclined to believe also the 3d quantum system to 
be integrable. If this is indeed the case is another a question worth studying. In this vein one 
could also try to generalize the Bethe ansatz for the XXZ spin chain to the constrained system 
of spin chains corresponding to the 3d case. Furthermore, we wonder whether the SUq(2) 
quantum symmetry of the XXZ spin chain carries over to the 3d melting quantum crystal. 
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A One-point function 

In this appendix we describe in detail the computation for the one-point function for the 
two-dimensional kink introduced in Sec. 14.3.21 

The starting point is the one-point function for the generating function of the kinks. 
Given the representation in terms of direct product over the spins (grand canonical ensem- 
ble) one finds: 

1 _ 7 n a x+l/2 



(Y(z)|^|Y(z)) 
Using the development 



-1+2E(-1)VJ*, 

k=0 



one can expand: 



(A.1) 



(A.2) 



1 oo N n K 

(Y(z)|^|Y(z)) = -i(- W ^ 2 ;, ? ) co + Y:E(-^ k ^ +m){N - k) l. 

N=0 ic=0 



-,k 2 /2 



Oft?), 



-w 



N 



(A3) 



The contribution of the N-particle kink is therefore 



m N (x) 



(Y N |T N ) 



N 



+ E 

k=0 



(A.4) 



where [^] is the ^-binomial coefficient 



(?;<?) 



N 



{q;q)x-k{q;q)k ' 



(A.5) 



Using the identity (see 11231 ) 



{q- N ;q)k 
{q>q)k 



(-l) k q Nk q-&, 



(A.6) 



one finds 



m N (x) 



-i N 



k=0 



(A.7) 



This represents a kink centered in x = N. For later convenience it is preferable to shift the 
center to x = 0: 



N 



m N (x) 



+ E(q- N ;q) k q m " kx 



(A.8) 



k=0 
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In the N — > oo limit one can use 

lim(q- N ;q) k q Nk = (-l) k q&, (A.9) 

N— >oo 

and find the exact expression for the shape of the kink 

m eo(*) = -. + E(-l)VV X - (A.10) 

To study the q — > 1 limit, introduce the function 

One can verify that satisfies the ^-difference equation 

Sf q {qQ = \-f q {Q, (A.12) 
and introducing the (^-difference operator as in [23 j: 

p ^O- /( fil{){° ' (A - 13 » 

one obtains 

(1 - (?) D,/,(0 = + ^r/,(0 • (A-14) 
In the q — > 1 _ limit, becomes the derivative with respect to £. Define 

/(x) = lim/ 9 (^ /lo sW). (A.15) 
Solving the differential equation for f(x), one obtains 

/(*) = A-^PlFT + ^^-^FT), (A.16) 



where E n (jc) is the exponential integral 



E n (x) = ^ — df. (A.17) 



Using the asymptotic expansion for n — > oo: 



t— — 7— / (A.18) 

x 1 + n/x 
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we obtain the asymptotic expression for f(x): 

/(*) = l^c ■ (A-19) 
It follows that the magnetization profile is 

m(x) = -\+f(x) = -^tanh(|) . (A.20) 

The shape of the partition is therefore given by 

u<p(x) = -2 J dxm^x) = -21og(2 cosh(f)) (A.21) 

or, in a more symmetric form (take u,p{x) =: y), 

e (x-y)/2 + e -(x+y)/2 = 1 (A22) 
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